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Abstract 

The  transition  probability  for  transfer  of   energy  from  one  optic 

to  two  or  three  acoustic  modes  of  vibration  has  been  calculated  for 

l2 
NaClj,   and  found  to  be  of  the  order  10       per  sec  at   iow  temperatures© 


1 


lo       IntrcxiuctioK 

In  a  crystalline  solidg   the  at  cans  oscillate  about  their  mean  positions  on 
the  regular  crystal  latticeo     At  not  too  high  temperatures^   these  oscillations 
may  be  considered  to  be  hamionieo     In  this  case^    if  we  Fourier  analyze  the  dis= 
placement  from,  equilibrium  of  all  the  atoros  in  a  given  crystal^,  we  can  describe 
the  3N  degrees  of  freedom  of  N  atoms  in  terms  of  3N  Foiurier  components  or 
"lattice  vibrations  "o     These  components  represent  the  thermal  motdon  of  the  sys-= 
tem  in  terms  of  vibrations  of  the  crystal  lattice  as  a  whoieo     (Strictly  speaking^ 
there  are  3N=5  such  modeSg   if  we  eiirainate  the  motion  of  the  center  of  mass  of 

the  system  as  a  wholep   and  any  rotation  about  seme  space  directiono)     In  practicej, 

fil 
following  Born  ana  von  Karman'='-c,   one  imposes  the   "cycxic  boundary  conditions'   to 


limit  the  number  of  acceptable  lattice  vibrations  to  3No     This  may  be  regarded 
W   as  a  device  for  obtaining  a  complete  set  of  lattiee  vibrations  in  terms  of 


[2]. 


which  to  describe  the  motiono  Use  of  the  cyclic  boundary  conditions  is  certain^ 
ly  permissible  if  N  is  sufficiently  largej,  and  if  we  study  phenomena  that  do  not 
depend  critically  on  the  surface  of  the  crystal o 

The  xattise  vibrations  in  a  polyatonie  crystal  are  of  two  typeSj,   namely 
acoustic  modesj,  whieh  correspond  to  translation  of  the  atoms  or  molecules  in  a 
gasj   and  optie  modes <,  which  correspond  to  intra:'nal  molecular  vibrations  and  ro<= 
tationo      (Thus  there  can  be  no  optie  modes  in  a  monatomic  crystal  like  Na  or  Neo) 
The  relation    between  freqiieney  (h/Zk  and  wave  momber  k/2m  for  these  two  types  of 
motion  is  shown  in  Figo  1  for  a  typical  line.ar   (one=dijnen^onal)  lattice-o     In  a 
three=  dimensional  lattice  there  will  be  several  optie  and  acoustic  branches,,   while 
the  limiting  wave  number  k^  will  in  general  be  different  for  different  directions 
in  the  crystalo 

The  existence  of  these  lattice  vibrations  as  stationary  states  of  the 
system  is  an  immediate  consequence  of  the  haormonic  character  of  the  oscillation 
of  the  individual  atoms  about  their  equilibrium  positionSo     In  fact  these  os= 
cillations  are  not  strictly  harmonic^  and  thus  the  description  of  the  thermal 


motion  of  the  crystal  in  terms  of  lattice  tibratloiig.  is  not  an  eatact  onej,   thoiigh 
it  is  a  very  satisf  aetory  description  wittiin  its  limitationso     If  we  expand  tha 
potential  energy  of  an  atom  in  powers  of  the  displacement  fifom  its  equilibrium 
position  in  the  lattice^   and  go  beyond  the  quadratic   (harmoniQ)  term  to  cubic 
and  perhaps  quartic  terms  in  the  displacement.,,  we  find  that  a  number  of  new 
physical  effects  can  be  explained^,   such  as  a  nori^^ero  thermal   expansions,   a 
finite  thermal  conductiidtyj   aad  the  po'ssibiiity  of  a  rise  in  the  specific 
heat  above  its  equipartition  values  at  sufficiently  high  temperafcureso 

If  we  take  into  accorunt  these  aiih^armoriic  terms  and  still  try  to  describe 
the  motion  of  the  atoms  of  the  erjstai  in  terms  of  lattice  '/ibrationSj,    these 
modes  are  now  no  longer  stationary  states  of  the  systemc,   since  the  anharmonie 
terms  in  the  potential  energy  induce  transitdons  beti^een  the  various  modes o 
Such  a  mechanism  is  invoked  in  the  e^lanation  of  the  thermal  conductivity  of 
ionic  crystals '"-J  ^  L  J  ^  f^r  exampieo     In  particular.,  by  reference  to  Plgo  1  we 
see  that  the  coupling  of  optic  and  acoustic:  modes  of  vibration  in  a  solid  is 
the  analog  of  the  excitation  and  de=>excitation  of  molecular  vibrations  by 
collisions  in  a  gaso     A  specific  application  of  this  process  is  to  the  caicu= 

nil 

lation  of  the  rates  of  various  chemic^'SL  reactions'--*  o  In  the  present  paper  we 
calculate  the  transition  probability  governing  the  flow  of  energy  between  optic 
and  acoustic  modes  of  vibration  in  a  crystal  as  a  starting  point  from  which  it 
might  ultimately  be  possible  to  attempt  to  caiciilate  the  rates  of  certain  chem= 
ical  reactions  in  liquids o 

For  simplicity  and  definiteness  we  shall  here  eonaidsr  the  case  where 

(1)  a>^  <     2c^  , 

so  that  an  optic  mode  can  split  up  into  two  acoustic  modes*     We  shall  .also  con^ 
sider  the  tr.ansition  into  three  acoustic  modeso     This  means  that  we  could  apply 
the  present  considerations  to  NaCl^   for  instance,,   and  we  do  in  fact  use  NcCl 
for  the  present  order=of=magnitude  caiculaticnso     However,,  it  must  be  stressed 
that  the  present  calculations  are  schematic,,    and  should  not  be  applied  too 


literally  to  any  given  substanceo 

2o   Harmonic  appro3djnation. 

We  consider  a  crystal  containing  2N  atoms<,  N  of  each  of  two  chemical  species 
A  and  Bp  in  a  regular  arrayo  The  displacements  of  these  atoms  frcan  their  equili= 
brium  positions  will  be  denoted  by  «  „  g  where  j  ■=>  Ac,  B  for  the  A  and  B  atoms  in 

a  given  basic  cell  which  is  labeled  by  n  ^  (n  si^^n  )o     There  are  thus  N  cells 
each  of  voltime  CX    g   the  n=th  sell  being  located  at  d  e  Now  we  expand  u,      in 
terms  of  a  set  of  normal  coordinates  q^,  ^  ^  <l^   <e(t)p  we  get 

»  1  kod 


Here  v„(kps)  is  a  unit  polarization  veetor^  s  ^  opto|,  aeoustoc.  corresponding  to 
the  respective  modes^  and  the  stanraation  over  the  wave  vectors  k  extends  over  all 
points  inside  the  first  Brillouin  zoneo  Note  that  our  normal  coordinate  q, 
has  been  defined  with  the  dimension  of  length,  not  (mass)    x   (length)|  in  this 

we  are  following  the  convention  in  this  field'-'^J^  L-^-J  o 

The  potential  energy  of  the  lattice  in  harmonic  approximation  is  given  by 

(a)  V     ^     7^^^  *  ¥^2) 

(3)  (b)  V^°^^     constant 


(c)        .     V^*^^^  {x/2!}  F".    o„  r°  „  A  ,      .„    „  u.     u 


where 


(U)  A.        .„    „s8Aoo„Cd[    .=  d„)     ^     A„,  „(d    „  =>  d    )  « 

since  the  potential  energy  depends  only  on  the  relative  orientation  of  the  cells 

r2i 

in  the  lattie@'--»o  _   ■ 

With  the  transformation  (2)|,  we  find  that 


(a)         V^2)         ^     IEv  ^  M^ 


l4 


where  aj(kss)  is  the   (angtaar)  frequency  of  the   (k^s)  mode  and  is  giTren  by  a  root 
of  the  secular  equation 

{=         ■=•  2  ik'od     I 

In  the  present   case  of  a  diatomic  lattice^   the  secular  equation  is  a  2  x  2 
determinant c    having  two  solutions  coCkgs)  for  eaeh  of  the  3N  Yaiues  of  k  lying  in 
the  first  Brillouin  zone|  these  two  solutions  correspond  to  s  "^  optop   acousto 
The  problem  of  solving  (6)  to  give  the  functional  form  of  «(kj,s)  and  of  the  fre= 
queney  spectrum  of  the  lattice  tabrations  was  first  discussed  by  Born  and  von 
Karraan'-  ^ §  for  any  specific  case  it  generally  ieadg  to  vary  lengthy  calculations o 
Here  we  shall  use  the  simplest  possible  approadmaticng   of  the   "Debj^e^Mnstein" 
form  (cfo  Figo   i)8 

(a)  Oritic  branch s         co  =*  oi       ^     independent  of  ^      1^1  £ 

(7;      ;        '  «  - 

(b)  Acoustic  branchg  ■(»  «  '2:.  jkfj,      |kj  <  lc.|     e.   ^  sound  velocityo 

Physically^   the  Einstein  form  (7^)  for  the  optie  braneh  corresponds  to  a  situ-- 
ation  where  the  internal  moleeuiar  vibration  is  only  weakly  coupled  with  the 
•external'   or  translational  motion©     It  is  a  good  appr-oxiraation  for  a  diatomie 
lattice  with  very  different  masses  of  ^he  two  ions|  in  this  casej,   howeverj, 
ca    »  coj.  -  c.k^o     The  Debye  form  (7b)  for  the  acoustic  branch  is  derived  on  the 
assumption  of  a  continuiiMr,   with  no  characteristic  atomic  structure o     We  see 
thus  that  the  present  form  of  the  relation  between  as  and  k  is  not  particular^ 
ly  adequate^  however^   it  is  certainxy  a  proper  form  to  take  for  an  exploratory 
calculation  like  the  present  one^  because  ariy  more  satisfactory  form  will  en= 


tail  much  more  labor  in  the  calailatione 


3o   Anharmonicitj 

If  we  take  third  and  f  oiirth  order  terms  in  the  potential  energy  into  account^ 
(3a)  must  be  replaced  by 


'^  (a)     V         «     V^°^  *  V^'^^  *  V^^-^  *  ¥^^^  *     •""> 


(8> 


(b)     v^^^  -  (L) 


BjH  Hn  n»n»«=J2^J''n«|,3»«n«»2.jfl  Hj9n«-j«»n«« 


rih) 


?1    \    ^=  ^" 


i^)     V  ^U.5''  ^J°j°°  j™°  ^=n  B'n"' n™"  "Jn^,j»n»^j"B"»  ^j»»"n«™ 


ijn  B,j«n«  H.joon»«Hj 


OlUin  OW 


If  we  transform  V^^^  and  Y^'^^  by  (2)s  we  find  that 


(a)     V^^^  ^ 


^ 


33N- 


=ka|k«s«|k»'  s»»  bCksik's^-k""  s»»  )  qj^^  q^,,g„  qj^„  g„ 


(b)     V^^'^  «   (rw)  >   ,      ,»    «,™     ™     ,  o«n    „,„e(kssk«s«^k»»  8»*|]£?«'B8«««)q>^  q. 


»s» 


k"  s  ™  "k'^s"""^ 


where 


( 


(a)     b(ks^°<">^ 


1 


"0    0   ^on     00  B  o      o„    „    c  „„    „„  V.  (ks)¥o  „(k"s")v.  „„  (k"" 


i  ^kod    4.  k»od   ,  *  k»«  »  d    „ 


(10) 


c(ks|°  °°) 


f  ^jn^^o^j 


.™    000  C-     „.„-,™    „™ir.(ks)'"'»v.„„(k'»«s™) 


exp 


V 


-l^'^n*  '""  ^  ^'^"in^n}- 


L2JsD] 


By  consideration  of  the  crystal  symmetry  it  can  be  shown.^ -^  ^  ^-^   that 
bCks^k's'ik""  s»"  )  and  c(ks^k»s«^k««  s"«  ^k«i'a")  differ  from  zero  only  if 


(11) 


(a)  k-s-k'+k""   s  0  or^a  vector  in  the  "reciprocal  lattice  t 

(b)  k  +  k"  *  k™  +  k«»  =3  0  or  Ko 


Further,,  h{kB^oco)s,   c(ks^ooo)   as  defined  in  (10)  are  independent,  of  N^  that  iSj 
of  the  total  size  of  the  crystaio 

It  is  rather  difficult  to  visualize  the  precise  me.aning  of   (lO)j,   and  so  we 
shall  present  a  simple  model  as  an  ixnastrationg  namely  a  monatomic  linear 
lattice  consisting  of  N  identisal  particles  spaced  on  a  line  with  mean  inter= 
atomic  separation  do     If  the  interaction  between  two  particles  distant  x  ^  nd  *  u 
is  9Ws>   then  if  we  expand  0(x)  in  a  Taylor  series  about  ndg  we  get 
^      (a)     V       -     V^°^  +  V^2)  ^  ^(3)  ^  ^(U)  ^   ,„„ 


n 


(12) 


i 


(b)  ¥''''^''s  constant 


,l^™.lF^°°(-d)(i=) 


'1  % 


^   n+m  n^         ^   n  n^m 


w  ^^'^-  I  E'-idli^"'"-^''!!)  I  (>.„™=V^  *(.-„-„„„. 


V 


(^'  v'^'-  ie:.i 


:::ir"(n.d)(ir)ru„«=^,)'' 


■  n  n=m 


We  wish  to  transform  m  to  normal  coordinates  q^j  then 


(13) 


VL_ 


1/2 


N 


^k  \  ^ 


■k 


(Note  that  v,  =  1  in  the  present  casso)  Let  us  suppose  that  each  particle  only 
interacts  with  its  nearest  neighbors.,  ioeoj,  m  -  N"  ^  lo  This  simplifies  the  ar= 
gumentj,  but  does  not  impose  any  restriction  on  ito  The  subie  term  is  now 


v(3) .  i 


m 


51"  L=n^l  l^kk^k"  Vk»\?'  '^k^t^" 


i(k*k»+k»«  )nd 


(lU) 


(e^^'^=l)(ei^"^^i)(e^'^'"'^=i)  *  (l=e-i^^)(l=e--^"'^)(i.e°^^""  ^)l 


Now  carry  out  the  sum  over  no  We  use  the  result 


[2].  [3] 


(15) 


=N      ifnd 


N  if     f d  ^  0     or     2nx  (an  integer) 

0  otherwise  „ 


which  is  precisely  the  origin  cf  condition  (il)o     In  this  way  wa  gefc 

/ 

(3)  ^    0M(d)     y=  rw    Y     ^(e"'^'^=iKe-^'^"'^-=i)f@'^"     -=l' 

2o3^jjV2  ^kk'k"  W«\«»<^®        ^^^®  ^''^®  -^^ 


(16)     ^ 


.    f.      ikd>..      ik»d^^,     iki^dvi   ^  „     „ 


3i 


1         

^72  Ckk»k"  ^(l^s^S^"  ^^k^k»^k 


00     s> 


where  the  only  contribution  comes  from  teiTTis  in  the  sura  satisfying  the  condition 


(l7a)  k  *  k"  •«■  k""    ^  0         or       2ioc  (an  integer)         ("reciprocal  lattice«)< 


Simiiarlys 
r 


(18)      < 


V 


c(k 


In  <^kk»k»'^  ^x'^s'^  3'^    5 


(k,k%k«sk»«)  qk\o\o,\™ 


o  o  o  o  < 


...)  .  tM^  ......  .,    J(e^'=^=i)...(.i'""<i=x) 


2  k  kB  k»»    k'"" 


^    /-,     "ikd.,,       (^     "ik^d-.! 
*•  (i=e        jo'oCl"©  )y    5 


where  the  only  contribution  to  the  sum  over  k"s  comes  from  terms  satisfying 


(i7b)  k  +  k'  +  k'"  -a-  k"     «     0  or       2iac  (an  integer)   ("reciprocal  lattic4)< 


=.  8  = 

It  is  Clear  that  the  present  argument  can  be  extended  to  any  order  of  neigh= 
hor   interactions  N'g  by  going  thr'ough  the  present  analysis  for  m  =  2<,  3<,  ooo 
guceessiveiyj  and  we  can  also  extend  it  to  thr'ee  dimensionso 

An  estimate  of  the  order  of  magnitude  of  b(ks|k»s«^k"'-'  s""  )  can  be  made  by 
■aaing  Kienens"  result"- -"'o  For  three  aeoustlc  modes  A,,  A%  A»»  he  gives  the  esti= 
mate 


(18)      bCAj-Ai^-A""  )  -  ayj'  M, 


<^\<^A.^^An 


'aeoust 


where  t  i§  the  Grueneisen  constant '-■^-'  defined  1^ 
U9]  Y    =     d  in  V 

which  equals  appj"oximateiy  2^   for  most  materiaiSo     The  result   (18)  disregards 
specific  crystal  structiire  and  polarization  effectSo     A  plausible  general issat ion 
of  this  to  the  present  case  where  we  want  b(0^Ai|A2)  would  be 

(20)  b(0^AijA2)     ^     (8/  yf  )    (M^p^M^,^^,/)"/^  T  <^^^^  /  e^  « 

It  is  even  more  difficult  to  give  an  estimate  of  the  quartic  potential  co'^ 
efficient 5   c  (0|;Al^ A2^ A3  )  o     On  purely  dimensional  considerations  we  shall 

^'^'^  b(6|AliA2r^  I  ' 

where  c  is  probably  in  the  range 
(22)  0,1    g     a    ^     10. 


lio       Phonon°phonon  coIllsioiMg  general  discugsion 

We  are  concerned  here  with  transitions  in  which  an  optic  phonon  Oj  of  frequency 
CO  and  wave  number  k^  goes  into  either  two  or  three  acoustic  phonons  Ai^  A2j,  A3  of 
frequency  coj^,^  at^^s  ^X3  ^^'^  ^^"^^  nraiber  k^^g  k^gs  1Ea3°  Thus^  we  wish  to  discuss  the 
two^phonon  process  lis 


(23)     Process  lis        0  =^  Al  *  A2 


p 


which  can  be  produced  by  the  cubic  potenti.al  V       j>   and  also  the  three=phonon  process 
Ills 

(2U)  Process  in§  0    •=>  Ai  *  A2  +  A3j, 

(u) 

which  can  be  produced  either  by  the  quartic  potential  V^  '^  or  by  successive  appli= 

cation  of  V^  o   "Graphs"  of  the  possible  mechanisms  of  these  transitions  are  shown 
in  Figo  2o 

In  studying  a  given  transition^  one  always  has  to  examine  the  conditions  of 
conservation  of  energy  and  momentumo  For  processes  IIj,  III  these  are 


(25) 


^o  "  ""Al  *  "A2 


%     "  "Al  *  "A2  *  "^AJ  ' 


^°      ^o  °  ^^M  *  ^A2^ 
(2o)    "x  y   ^  0  or  K,at vector  in  the  re° 

■^°     ^o  °   ^"Al  *  ^2  *  =A3^  ciprocai  lattices 

Notice  that  on  account  of  the  periodic  structure  of  the  lattice^,  momentum  has  only 
to  be  conserved  modulo  a  vector  K  of  the  reciprocal  lattice^  processes  with  K  f^  0 
are  caxled  "IMklapp"  processes '^ -^  ^ '•-' 0 

In  calculating  the  transition  probability  for  a  given  process^  we  use  the 
quantum  mechanical  rather  than  a  classical  descriptionp  and  in  particular^  we  use 
a  field  theoretical  description^ -"^ '•J  ^  *= -'  which  is  convenient  and  well  knowne 


10 


Thus  we  need  the  transition  matrix  elements  for  these  process? 
The  transition  matrix  element  for  process  II  is 


(27)     VV,    ^.      (F|V^3)jo).    (V2)3/2b(Mi|A2l 


II 


L%t^acoust%'"Al"A2J 


where  n  ^  n.sp  n.^  are  the  roumbers  of  phonons  of  type  Oj  Al^  A2  before  the  colli = 
siono  Similar ly^  the  transition  matrix  element  for  process  III  is 

r(3)(TVThT(3)| 


III 


(28) 


,.,.,2     ^(0|Al|A2iA3)    I      %faAi^^)UA2^i-)faA3^-^>^- 

J 


hi  N 


-  opt.  acoiist  o  Al  A2  A3 


c(0|Al|A2^A3)  -  c(0|Al|A2.^A3}   * 


2b(0_8Ai^3A)b(A|A2^A3)(n.*l) 


\coust^%°'"Ar'*i  ^  ^r 


Recall  that  whixe  monentimi  is  conserved  in  both  steps  of  the  transition  IIIBj, 
namely 


(29) 


Al  *  A 


A  ~>  A2  ■*■  A3c 


yet  energy  is  not  conserved  in  the  indivldiii.ai  steps  of   (29)j)   though  it  is  of  coiarsa 
conserved  overall o     Thus  the  energy  denominators  in  (28)  do  not  vanish o 

In  order  to  find  the  transition  probabiiityg  we  still  have  t*D  find  the  den^ 

sity  of  final  states  per'  unit  energy  range  o     For  pro<e.ess  II  this  is 


(30) 


dp(ll)/dE 


"^^/^M    Na  d  cose^^d0^^k2j   di^^ 


(2n)^dE 


(2j!)- 


dE 


dk^^/dE     -     1  /J   Kc_p(l=  cos©^.,)V    s 


where  ©^, „  Is  the  angle  between  the  vectors  k..^  and  k.«c,    0.^^,  is  the  vo±tiine  per  imit 
cellg    and  c.   the  velocity  of  sound  of  equation   (7b) o     For  process  IHj,   if  we  treat 
^Al^  ®A?^  ^A2  ^^  independent  variablesg  we  hav© 
/ 


(31) 


dp(in)  Nn,d\i      Na,^k22dcos9^2d0^2       dk^ 


< 


dE 


A  2 


dE 


{2ny  {2%y 


dE 


where  G„-.  is  the  angle  between  k.«  and  k.^  o 

Now  the  transition  probability  per  tmit  time  for  a  given  proces 
III  is  given  by  the  formula 


s  A 


II. 


(J2) 


P^       ^     (2nA) 


V« 


dpCA)/dE|         A-     11,    Ills 


where  )    „  denotes  summation  over  an  final  states  consistent  with  conservation 

of  energy  and  momentumj,   that  is,,   integration  ovear-  d  cos9,^   d0.„  in  case  11^   and 

■a  ,     . 

0¥er  i-^k.-^d  cosS,^  ^^A?  ^  sase  IIIo     We  note  that  P  is  a  well=defined  physical 

quantity  in  that  it  is  independent  of  \  ioeot,   of  the  size  of  the  crystaio     It 

doesj   however^   still  depend  on  the  magnitude  and  directiai  of  the  incident  m©men= 

turn  k  e     This  is  partly  a  genuine  effect g  but   to  some  extent  due  to  the  Einstein 

foriTi  of  the  frequency  spectrum  (equation   (/a))o     At  a  later  stage  we  shall  average 

.■file  transition  probability  over  ail  possible  incd.dent  momenta  and  also  over  ail 

possibilities  of  linklappo 


We  wish  to  make  a  detailed  calculation  of  the  transition  probability  for  pro= 
cess  II  of  equation  (23 )o     Let  us  begin  by  considering  the  restrictions  imposed 
by  conservation  of  energy  and  m^nentumo     Equations   (25-11)  and  (26=11)  give  four 
conditions  for  the  six  unknown  con^jonents  of  k.^  and  k.^   so  that  there  are  two 
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imdelentiined  parameters  (whieh  were  taken  as  t.he  two  angxes  ®,js  Pa-j  i"  "the  dig- 
oission  of  equation  (30) )«     Cieariy  k^„j,   kj^-  define  a  plane  =  say  the   (xy)=plane| 
k    lies  in  that   planec,   and  may  be  taken  to  define   the  x-aadso     Thios  the  ajdjtmth 
giving  the  orientation  of  the  (xy)=pian8  in  space  is   one  of  our  two  undetermined 
parameterso     The  second  parameter  is  less   tri^iai  and  more  specific  to  the  modelg 
so  let  us  now  use  the  explicit  form  (7)  of  the  relation  between  to  and  ko     If  we 
write:  q     for  k  o   or  for    (k=  K),,   depending  on  whether  or  not  there  is  an  IMciappj 
then  conservation  cf  momentum  and  energy  impose  the  following  requirements  on 

Mai  ^^'^  -A2° 

^33)  k^i  *  k^2    ^     2o  ^  ^Al  *  ^A2  ^^  ^V^A>° 

Thusc,  the  endpoints  of  k.,  and  k.^  lie  on  an  ellipse  (strictly  speaking^,  on 
an  ellipsoid  of  revolution)  o  The  situation  is  shown  in  Figo  3s  and  we  see  that 
our  second  undetermined  parameter  may  conveniently  be  chosen  as  n  ^  cos  a,,  rather 
than  as  cos  ©.  ©  The  ellipse  is  characterized  by  its  eccentricity 

(3Ua)         e  -  OB/OX  -  (%V%^ 

and  its  semi=-major  sodB 

(3iib)  a     s     OX     «     %/2e^  e 

In  terms  of  these  we  can  readily  obtain  the  results 

/ 

k-,     ^    a(l  + 


H 


i=cos  e,„    -   ^^K®J    « 


^  l=e2^.2 


\ 


We  have  to  investigate  the  possibility  of  Umkiapp  processeso     This  is  dis= 
played  graphically  in  Figo  h  for  typical  values  of  the  magnitudes  of  k  g   k.=,fl    and 
k.po     We  see  that  in  addition  to  the   onor^^^'  '  process^   there  is  just  one  possible 
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type  of  Iteklapp  in  the  present  sitnatioc,* 

Let  us  now  eonibine  eq'uatiorB(30)j,    (32)s  and  write  the  transition  probability 
'.or  process  II  as 

(3&)  bfyj/djii,     ^  i»    0     ^     ■■  P\>L  i=rosS  '  ^ 


^11 '5 


Where  (%j^  and  hence  ?» ,  is  regarded  as  a  function  of  p.  through  C?)^    (27)^  aM   (35)< 

We  are  interested  in  the  total  transitioai  probability 
m 

^Il^io^     ^     J 


^jj/dp.    dju  o 


Note  that  in  general  the  limits  of  integration  will  not  be  =1  <  ji  <  lo     We  know 
that  k .   5   k.^  <  kjp  which  gives  the  Jlmits  =m  <  t^  <  %  where 

(  \ 


(38)  in  *-  Min'ds    (i/e)   (2©^,/©  =i)V    «    Mn^lj    CV%^^^^©^^^^    " 

I  J  [  J 

Combining  C20)<,    (27)5    (30)  and  (37)s  we  obtsdn  the  estpressim 

ti  n,,TV jB^^^l)  (B^g-M)      f"  4^co^g         acos9^3L     ^ 

We  must  note  that  P,j  depoads  on  the   order  of  Ifekiapp  through  q^y  which  is 
equal  to  k^  ©r  k    =  K»   depending  on  whether  we  are  dealing  with  a  normal  process 


(N)^  or  an  Umklapp  process  (U)c,   "whieh  eerresponds  to  the  case  n    ^  1  of  I^lgo  iio 
If  we  write  k     =  ^Kis  then  ow  approadmate  finax  esqression  fcr  ^ttC^q)  is 

'PjlCq^')  -  Fn(K.A)  "  P^TL  (K)n^(n^*l)(n^2*l)|  A-     N,U 


(uol 


po     .     }%^o^^A^^^ 


ii32n(M     ^M:         .)^^^  .    , 

^  opt  aoousfc 


\ 


(UO  eont'd] 


AK)    -  iK^/2=K)(x=  Kf/6^  y[ 


<  Yk  <    I* 


In  practicej,   we  are  most   interested  in  the  total  transition  probabiiityv,   repre«= 
senting  the  aTsr-age  of  P^.^(k  „A)  owet  an  possible  values  of  k    g 

f  J      f 


V 


J 


The  actual  ntrniericai  values  used  for  the  various  parameters  are  listed  in  Table 
Ij   and  gi-ve  the  final  result 


(a2) 


II„av^   ^ 


< 


-     5o8  X  10      per  sec 


jn^'ZiZ,^.,         -    3-2  -  lO^Ser  =e, 


'o'^'mi^'^"A2^ 


for  the  mean  transition  probability  per  Tjnit  time  of  one  optic  phonon  splitting 
up  into  two  acoustis  phonoas  under  process  H  of  equation   (23) » 


5e       The  three°phonon  process  III 

The  purpose  of  this  section  is  to  sketch  the  caictilation  of  the  transition 
probability  for  process  HI  of  equation  (2U)  and  Figo   29     The  caiculaticn  will 
not  be  carried  out  in  the  same  detail  as  that  of  the  preijious  section  for  the 
process  IIj,   our     ofo^'eet  being  to  indicate  that  process  III  merely  furnishes  a 
small  correction  to  picoess  IIo     In  particular^  we  shall  not  go  into  detail  in 
examining  the  dependence  of  the  matrix  elements  on  k     ^   k.^p  ^AP^  =A'(S  because 
our  limited  knowledge  of  the  coefficeients  b(0|Al|A2)  and  e(O^Ai|A2|A3)  does  not 
warrant  it« 

From  (27 )s    (31)  and   (32 l^   the  transition  probability  for  process  III  is 


6 


Vfe  now  have  to  integrate  over  d'k.^  d  cos9,2d0,^o     Let  us   recall  that    |c(0|Ai|A2|A3) | " 
depends  on  these  quantities  in  some  unknovm  but  complieated  way  (cfe   Section  3)s   and 
hence  let  us  replace  it  by  an  average  vaiue<,     Similarly^    let  us  put 


(Uii) 


J    (2n)\^    J     ^^=^°^^23^  '^l     ' 


which  gives 

f  .s       .       f      ^       ^^•^c!^(°^A^M2gA3)ilf^nJn^^.i)(n^g.l)(n^y.i) 

9216  H  M^p^M^^^^^^^  e^  cc^  o^ 

This  may  be  compared  with  the  corresponding  result   (UO)  for  P,^(q  )o     At 
sufficiently  low  temperatures  we  may  put 

(UO)  n^^     -    n^2     "    ^3     "     °^ 

as  relatively  few  acoustic  phonons  will  be  excited  irdtiaiiyo     Then,,   disregarding 
numerical  factors  of  order  one,,   we  have 

U.^  ^in^io^  |e(0|Ax|A2|A3)|J^  ^ 

Khl)  fc     t„    \        ^ 


^V^  |b(0|Ai|A2)|?  -^^^coust^M 


av 


If  we  now  make  the  substitutions  (20) ^  (28)  for  bCooo)^  ^{°'">)s   and  put 

nr+  1  ■, 


^'^o-^Ar'^r^  "^T.  '\coust  "^A  ^cotist 


then 


^mW'  (i^a/iOO)^ 


^^9)  =rnin  --  ^"13       ^      Ooi^^  a  ^  id 

this  ratio  is  certainly  very  mueh  less  than  onso 


7o   Dissussioa 

The  central  numerdcaL  resiilt.  of  the  present,  cajiculaticn  is  equaticn    (^2)3 
which  give§ 

(50)  P^-        <-  iO^^  see°^ 

at,  sufficient.iy  iow  temperatarsa;   (T  «  t^Kiiw/ic,   so  that  nu^  ^  n^^  ^^  0  in  the  mean)o 
We  may  con^are  this  with  oowj  whish  gives  aome  measixre  of  the  mean  "soiiision  fre^ 
quency"  of  atoms  in  the  crystalo     On  this  basis,,   therefore^   we  may  say  that,   the 
efficiency  of  thermsii  (aeoustis)   vibrations  in  de-exciting  an  optic  mode  is  of 
order 

(51)  ^n,aT^     ^    2  X  xO"2  . 

Simixar  mimerical  vaxnes  have  been  obtained  eariia'g   both  caassicaxiy  and  quantum 
mechanicailyn   by  Pauxic,   Born,   Biacicman,,   Barnes,   Brat.tain  .and  Seitz^o   and  Huang*-  -' 
in  calculating  the  widths  of  infrared  dispersion  lines  in  crystalso     The  result 
(5i)  may  axso  be  eompared  with  the  same  general  range  (10      =  lO     )  fcr  the  effi= 
ciency  of  de^excitation  of  moxeeaiar  vibrations  by  collisions  in  gases '^"'o 

A  number  of  assumptions  have  besn  made  in  the  caxcuxationo     One  of  these  is 
that  the  perturbation  expansion   (6)  conva"geSo     This  seems  to  be  entirexy  satis- 
factory! one  knows  that  the  cubis  and  higher  ordesr-  terms  in  the  pot-ential   are  not 
important  at  reasonably  1.0^  teiper'atureaip   and  this  is  certainly  borne  oiat  by  the 
calculations  of  Section  6© 

The  present  choice  of  frequency  spectrum  is  reasonable  far'  an  exploratory 
calculation  like  the  present  one5   bat  it  is  certainly  rot  particularly  satis= 
factory  for  a  substance  sueh  as  NaCio     It.  would  be  adequate  for  a  lattice  of  two 
kinds  of  atoms  having  very  different  masses<,   or  for  molecular  crystals    (which 
would  be  the  case  of  real  interest  here)o     In  thos©  caseSg   however,,    one  would 
generally  have 


(52)  ^0^%.  '^  ^s  s^  »  2s 

(n*l) 
so  that  n=th  order  transitionSj  and  terms  up  to  V^  in  the  potential   energy 

expansion  {Sa)^   would  have  to  be  considered  in  the  calculationo     This  would  give 

rise  to  technical  problems  in  caxcuiating  transition  probabiiiliesj,  but  more 

particularly  to  considerable  new  diffisiilties  in  estdmating  the  magnitude  of 

these  terms  in  the  potentiaio 

Regarding  the  detailed  expression  used  for  the  coefficients  of  the  cubic 
and  quartic  potentials^   equation  (20)  and  especially  (2l)'=(22)  can  hardly  be 
regarded  as  a  satisfactory  approximaticBi  in  the  present  caseo     More  detailed 
calculations  of  the  potential  certainly  appear  to  be  called  foro 

One  may  ask  why  the  present  calculation  has  been  referred  to  a  very  approxi^ 
mate  three-dimensional  model  rather  than  to  a  better t,   one-dimensionai  oneo     The 
reason  for  this  is  that  in  three  dimensions  there  are  relatively  many  more  high 
frequency  modes  than  in  a  one-dimensionai  lattice^  (N(cD)diio  f^  dm  In  one  dimension^ 
CO  dw  in  three  dimensions )5   and  this  is  probably  more  important  than  the  detail  of 
the  frequency  spectrumo 

I  should  like  to  thank  Dro   Ta=You  Wu  ar^d  Dto   G«  R»  Allcock  for  discussiaiSo 
'?stt  of  this  work  was  carried  out  at  the  National  Resesrch  Cofjncil^   Ottawaj, 
Canada,     I  should  aiso  like  to  tlmnk  DTo   Go   Herzbergc,  Director  of  the  Di -Bis ion 
of  P\ire  Physicsj  for  his  hospitaiityo 


Table  lo^    Numerical  values  of  various  eCTist^ants  fog"  NaCl 

Debye  temperatiireo        ®q  281  \ 

(n^  ^  k®j^/h  3o67  X  10^^  sec°^ 

a>  7ol6  X  10  -^  sec 
o 

interatomic  distance^   d  5»t>3  x  10°     cm 

Grueneisen  constantj,  y  1«61 

Debye  mean  soimd  velocitjc,  c.                                 2«73  x  10^  cm/sec 

\coust  9o71  X  10=23  gn, 

n^  io?85  X  10°^^  OT^ 
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optic  mode   ^ 


"^'^ont  ~  internal  fre- 
quency of  vibration 


I  /acoustic  mode<^ 


3tj     ,  =  0,  as 
^     acoust    ' 

translation  cor- 
responds to  zero 
restoring  force 


Solid 
(Debye-Einstein  model) 


Solid   (exact) 


Gas 


Figtire  1.     Modes  of  vibration  of  a  one-dimensional  crystal  lattice 


Process  II 
0-»>Al  +  A2 


Process  IIIA  Process  IIIB 

0  -»  Al  +  A2  +  A3 


Figure  2.  Phonon-phonon  collisions  (time  is  measured  increasing  upwards) 
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(acosa,  bsina) 


Figure  3.  Energy-momentum  conservation 
in  two-phonon  collisions 
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(i)  Normal:    no  Umklapp 
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'li)       n    =   1  :   Umklapp  is  possible 
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(iii)     n    =  -1   :   no  Umklapp  possible 


<— • 
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I 

2k 
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(iv)        n     =  tl:    no  Umklapp  possible 


'  Figure  U.     Umklapp  possibilities  in  two-phonon  collisions 
We  take  Ik^J    ,    1^2 1  ~    ^-^  k^^     '        ''^ol  ^  ^'^  V 
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